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Abstract Two and a half years prior to China's M7.9 Wenchuan earthquake of May 

2008, at least 300 million metric tons of water accumulated with additional seasonal 

water level changes in the Minjiang River Valley at the eastern margin of the Longmen 

Shan. This article shows that static surface loading in the Zipingpu water reservoir in- 

OO duced Coulomb failure stresses on the nearby Beichuan thrust fault system at <17km 

Cn depth. Triggering stresses exceeded levels of daily lunar and solar tides and perturbed 

^_^ a fault area measuring 416±96km . These stress perturbations, in turn, likely ad- 

r^ vanced the clock of the mainshock and directed the initial rupture propagation upward 

^ I towards the reservoir on the " Coulomb- like" Beichuan fault with rate-and-state de- 

JL pendent frictional behavior. Static triggering perturbations produced up to 60 years 

flj (0.6%) of equivalent tectonic loading, and show strong correlations to the coseismic slip. 

(3JQ Moreover, correlations between clock advancement and coseismic slip, observed during 

jyT the mainshock beneath the reservoir, are strongest for a longer seismic cycle (lOkyr) 

i^ of M>7 earthquakes. Finally, the daily event rate of the micro-seismicity (M> 0.5) 

^ correlates well with the static stress perturbations, indicating destabilization. 
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1 Introduction 



> 

in 

iy-\ The Wenchuan earthquake of May 12, 2008 occurred in the Sichuan province of the 

y—{ People's Republic of China. This M7.9 event ruptured along the border of the Longmen 

O^ Shan margin of the Tibetan plateau in the West and the Sichuan basin in the East. 

r — , The earthquake's nucleation point was at Long=103.364 and Lat=30.986 at about 16 

(*^ km depth. The epicentral error is 5 km and the focal depth estimation error is 10 to 



O 15 km (Huang et al. 2008 1 



The Beichuan fault system, consisting of listric and NNW dipping reverse faults. 



broke 250-300 km parallel along the Longmen Shan thrust belt (Fig. 1) ( 
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et al. 19951. This intra-continental region has been extensively studied both prior to 



and after the M7.9 earthquake in 2008, including a) paleo-seismicity studies (Burchfiel 



et al. 1995| [Densmore et al. 2007) [Zhang et al. 2009[ [Zhou et al. 2007) [Burchfiel et 



al. 1995 



Chen et al. 19941 



b) instrumental recordings of the seismicity prior to and 
after the mainshock (1970-2009) ( ,Hu 2007 Lei et al. 20081, c) inversion data analyses 
of teleseismic body waves of the mainshock ( Wang et al. 2008a| Zhang et al. 2009] ), d) 
studies on coseismic ground deformations { Lin et al. 20081 Zhang et al. 2009). Finally, 



the trigger mechanism of the 2008 M7.9 Wenchuan earthquake has been debated since 



the earthquake's occurrence (Klose 2008 Lei et al. 2008). Some studies suggest that the 



mainshock might have been triggered by pore pressure diffusion within the earth's crust 



resulting from a nearby artificial lake, the Zipingpu water reservoir (Lei et al. 2008 



Ge et al. 20091. Other studies reject the hypothesis of triggering due to pore pressure 
diffusion ("Deng et al. 2010 Gahalaut and Gahalaut 20101. It can be anticipated that 



the earthquake cycle was already in its late stage in 2005 and close to failure conditions, 
because the Wenchuan earthquake ruptured in 2008. But, did the surface loading affect 
the earthquake cycle and the initial rupture? How many years of equivalent tectonic 
loading would the artificial loading advance the clock of the mainshock and how much 
of the coseismic slip would it produce? 

This study shows that observations and data modeling support the initial argument 
that lithostatic stress changes and the poroelastic response of the earth's crust due to 
the weight of the Zipingpu reservoir on the earth's crust most likely triggered the 2008 



M7.9 Wenchuan earthquake (Klose 2008 1. The article provides further evidence that a) 



surface loading due to water mass accumulations within the Zipingpu water reservoir 
in the Minjiang River Valley between 2005 and 2008 induced Coulomb failure stress 
changes in the earth's crust and b) triggering stress perturbations beneath the artificial 
lake likely advanced the clock of the Wenchuan earthquake, while affecting the initial 
rupture propagation. Furthermore, the water reservoir impounding generated static 
triggering stress perturbations in the earth's crust beginning in 2005 that biased daily 
stress alterations due to tidal elongations of the moon and the sun. Furthermore, it 
could have changed the natural earthquake cycle and advanced the clock of the 2008 
Wenchuan earthquake by several decades. 



2 Data 



2.1 Seismicity prior to the mainshock 



Instrumentally recorded seismicity with magnitudes up to 5.0 was observed in this 



region before 2004 (Yang et al. 2005 Liu 2007 Hu 2007 Yao et al. 2008 Huang et 



al. 2008[ ). Paleoseismic studies also show evidence of Quaternary reverse faulting in the 
Longmen Shan region. Radiocarbon analyses ( C), for example, indicate that the last 
major earthquake (M<7.9) might have occurred between 4 and lOkyr ago ( [Densmore et 
al. 2007 Zhou et al. 20071. Moreover, deformation measurements suggest the Longmen 



Shan is a transition part between both a stable continental region (Sichuan Basin) with 



strain rates < 10" 
rates > 10~^yr^^ 



-10, 



yr 



and an active continental region (Tibet plateau) with strain 
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Fig. 1. Geology and tectonic situation in the Longnien Slian region (a) the Beichuan fault 
system that ruptured during the May 12 2008 M7.9 earthquake with a coseismic slip 
s ( |Wang et al. 2008a[ l (b) and likely due to static Coulomb failure stress perturbations 
ACFS (c) as a result of the surface loading in the Zipingpu water reservoir (see Fig. 2). 
The principal lithostatic stress elipsods show the orientation of reverse fault contitions. 
The epicenter and the epicentral error of the mainshock are indicated by red circles. The 
geology is as follows: Qt-Quaternary, Cr-Cretaceous, Jr- Jurassic, Tr-Triassic, Pal-Paleozoic, 
Prot-Proterosoic, ProtG-Proterosoic granite, BM-Boashan Massif, PT-Pengguan Massif, 
WMT-Maowen thrust fault (WMT) and BTl/3-Beichuan thrust fault | |Burchfiel et al. 1995| 
|Jia et al. 2006| l. 



The area where the mainshock nucleated tends to have low horizontal and vertical 
tectonic deformation rates of <1.0±1.0 mm yr~ , resulting in less seismic activities 
when compared to a) the seismic hazard regions in the Quiangtan block (Southwest) 



and the Kulun Mountains (Northwest) and b) the Mesozoic deformation period ( Burch 



fiel et al. 19951 [Densmore et al. 20071 [Zhang et al. 2009[ [Zhou et al. 2007] [Burchfiel et 



al. 1995| Chen et al. 19941. High shear- wave velocity structures and Bouguer-gravity 
structures indicate a mechanically strong mid crust (10-20 km) and a cratonic-like litho- 
sphere of the Sichuan Basin (Fig. la), in particular, in the Yangtze block along the SW 



Longmen Shan at the border to the Songpan-Garze block ( Yao et al. 2008 Burchfiel et 



jal. 19951 ). ^^o ^^ ^^- suggested a more unstable lower part of the crust (>20km) along 
the Longmen Shan which could be explained by low shear-wave velocity structures 



(Yaoet al. 20081 



2.2 Seismicity associated with the Zipingpu reservoir impoundment 

Between October 2005 and May 2008 shallow seismicity patterns (<10km) were ob- 
served on the Beichuan fault system along a <20 km long part of the lower Minjiang 



River Valley parallel to the extension of the Zipingpu water reservoir ( Lei et al. 2008 



Klose 2009 1. Such shallow seismicity could be indicative for earthquakes in a stable con- 



tinental crust ( [Klose and Seeber 2007| ). However, Lei et al., reported that the seismicity 
change might be Eissociated with a continuous water accumulation by the Zipingpu dam 
in the Minjiang River Valley almost three years prior to the rupture of the mainshock 



(Lei et al. 2008 Klose 20081. At the northern end of the valley, the Minjiang River 



flows into the Chengdu plain of the Sichuan basin. The artificial water reservoir ex- 
tended sub-parallel to the Beichuan fault segment BTl about 1.5 km westward and 
BT2 less than 1.0 km eastward (Figs. 2 and 1). After continuous impounding, the wa- 
ter level peaked twice in October 2006 and October 2007 at water reservoir capacity 
of 1.10 10 m (upper water level) which is equivalent to a mass of 1.10 billion metric 
tonnes (Gt). Between October 2007 and May 12 2008 the reservoir was, again, season- 
ally drained to a remaining volume of 0.32 10 m (lower water level). Lei et al. first 
mentioned that the seismic events of magnitudes ranging between 0.5 and 3.9 might 
have illuminated a destabilization process on BT during the loading period of the wa- 
ter reservoir. Prior to the surface loading, the correlation coefficient between the water 
level and the daily event rate (383 days) is -0.18 (Fig. 3). Since May 2005 - the start 
of the flooding season - a positive correlation exists between the water level (mass of 
water) and the micro-seismicity beneath the water reservoir. With a coefficient of 0.69, 
the correlation is strongest between 2005 and 2006 during the first 366 days of loading. 
Between 2006 and 2007 the correlation is weaker with a coefficient of 0.62. Thus, it 
could be anticipated that the earthquake nucleation processes started already in 2006. 



2.3 Coseismic slip of the mainshock 

Numerical inversion data of teleseismic body waves (P-waves) show that the mainshock 
described a complex rupture process on the Beichuan fault segment BTl (Fig. lb) 
(jZhang et~al. 2009: Wang et al. 2008a[). During the first 15-20 seconds, a reverse fault 



focal mechanism dominated at and above the nucleation point of the mainshock ( Zhang 



et al. 2009 Wang et al. 2008a I. The rupture propagated with a coseismic slip of up to 



4 m on BTl upwards toward the Minjiang River Valley. Thus, it propagated directly to 
the surface loading area of the Zipingpu water reservoir. After 20 seconds, the rupture 
propagation process changed to a right-lateral NNE strike-slip mechanism NE of the 



Minjiang River Valley ( Zhang et al. 2009 Wang et al. 2008a I 



3 Methods 



3.1 Triggering stress perturbations caused by surface loading 



An exact first order solution of stress states below the surface loading area (lower and 
upper water reservoir) in the Minjiang River Valley (Fig. 2) is according to Boussi- 



nesq's classical solution (Boussinesq 1885 Love 19441. Boussinesq's solution for point 



loads is based on the assumption that the modulus of elasticity is constant within a 
homogeneous 3-dimensional half-space (earth's crust). Moreover, the principle of linear 
superposition is also assumed to be valid. For the given problem, surface areas of the 
lower and upper water level with arbitrary geometry were discretized in A'' Oj x 6^ m 
elements (e.g., 50x50 m ). With respect to the surface elements, A. Love developed 
a method to analytically determine stress states at a depth z beneath any uniform 
un/load Li with the area A, = a, x hi (Love 19441. This analytical solution is free of 
any elastic constant and can be defined for each 2-D surface element i = 1, 2, 3, . . . , A: 



o"L,i(-^"ij, z, ai, hi) = Li(Ami, Ai) 



f{z,ai,bi) 
2-K 



with 



(1) 



f{z,ai,bi 



arctan 



R, 



+ 



Ri 



+ 



bl + : 



(2) 



and where Ami is the mass change and Ri = \/{ai +bi + z"). f{z, ai, hi) describes 
the fraction of Lj in the vicinity beneath Ai and ranges between and 1. The un/loads 
Li induce a positive /negative vertical stress alteration Aa^ = Aaj^ superimposed over 
all (J[^i, whereas cr < is compression. L also changes the horizontal stress components, 
due to the elastic response of the crust (Hook's Law): Aa^ = Aai^ and t^/l-i/ Aai^ < 



^<^\,2 < '^'^L (McGarr 19881, where v is the Poisson's ratio. It is likely that Aa\^2 > 
vji-u Aaj^, because horizontal principal strains t\^i are very small and can be assumed 
to be ei 2 = in 11-19 km depth, where faults are generally locked in the Longmen 



Shan region ( |Wang et al. 2008b I. 

Table IT] shows the 3-diniensional model quantities that were taken into account to 
determine the lithostatic stress perturbations in the earth's crust. Stress states were 
determined on 7x4 km large fault elements of the double-listric Beichuan fault system 
(BTl, 2, and 3) near the artificial reservoir in the Minjiang River Valley. Thus, (Jj^^i 
and the resulting static Coulomb failure stress ACFS: 



ACFS = Art + tan (j) {Aa'n 



(3) 



were calculated on each fault element with respect to the lower and upper wa- 
ter level in the Minjiang River Valley (Fig. 2), whereby a'n is the effective normal 
stress. ZiC-FS-values were used for further estimations of the clock advancements At 
with respect to 4 kyr, 7 kyr, and 10 kyr earthquake cycles of the fault system (this 
methodology is described the next section). 



Pore-pressure diffusion most likely played a major role for the micro-seismicity 
(M<4) that was observed beneath the water reservoir, according to the study of Lei et 
al. 2008. Thus, drained conditions might have existed in shallower depth (< 10km) in the 
Paleozoic rocks and the permeable Mesozoic rocks. This suggests that the poroelastic 
response due to the lithostatic stress alterations of the surface loading played a minor 
role under these geological conditions. Thus, the Skempton coefficient B was assumed 



to lie between 0.5 and 0.7 (Terzaghi 1938 Biot 1941 Rice and Cleary 19761 



Ap = -BAa, 



(4) 



where Aa is the change of the mean stress. 

In undrained conditions (>10km) the pore pressure p increased due to the poroelas- 
tic effects and destabilized BTl and BT2 within the low permeable Proterozoic granitic 
rocks. The poroelastic response tends to increase with both depth (confining pressure) 
and fracture density ( [Lockner and Beeler 2003b| . The Skempton coefficient B, which 
describes the strength of the poroelastic response, was assumed to range between 0.7 
and 0.8 in 10-23 km depth. It is also assumed that B is additionally amplified up to 
0.85, due to a higher fracture density (Lockner and Beeler 2003b) near the intersection 
regions of BTl and BT2 in 15-17 km and at the South end of BTl in direction to the 
intersection with the Wenchuan Maowen thrust fault WMT (Fig. 1). Furthermore, the 
Mohr-Coulomb failure criteria indicates that the poroelastic response dominated only 
on steep dipping fault segments (>60°). Geomechanical parameters, which were used 
to determine triggering stress perturbations, are summarized in Table l!] 



Table 1 Model quantities for determining the Coulomb failure stress on the Beichuan fault 
system due to the surface loading in the Minjiang River Valley (Fig. 2). 



Gcomechnaical quantity 



MeanitSME 



Characteristic element on the fault 



Poisson's ratio u: 0.25ib0.05 

Skempton's coefficient B (drained): 0.60ib0.10 

Skempton's coefficient B (undrained): 0.75ib0.05 

Skempton's coefficient B (undrained): 0.85it0.25 

Rock friction angle </>: 28±2° 

Rock friction angle </): 28it3° 

Rock cohesion cq: lOitlO MPa 

Rock density p: 2700±50 kg m'^ 

Young modulus E: 75±25 GPa 

Shear modulus G: 30±10 GPa 

Lithostatic stress regime: cti > (T2 > 0-3 

Horizontal lithostatic stresses: ci and 1T2 

Gravitational lithostatic stress 0-3: pg^ 

Fault dip angle 9: 20-60±l° 



Mesozoic rocks (<8.5km depths) 
Mesozoic rocks (<8.5km depths) 
Proterozoic rocks (>8.5km depths) 
faults zones (>8.5km depths) 
Mesozoic rocks (<8.5km depths) 
Proterozoic rocks (>8.5km depths) 



reverse fault regime 

vertical stress 
BTl, BT2, and BT3 



3.2 Clock advancement of the mainshock caused by surface loading 



Paleoseismic analyses show evidence of Quaternary reverse faulting in the Longmen 
Shan region which indicate a M7 (or M8) earthquake recurrence interval of 7±3 kyr. 



( Densmore et al. 2007 Zhou et al. 2007 1 given horizontal and vertical deformation rates 



of e = l.Oil.O mm yr (5±5 10 yr ). The rate-and-state dependent friction law 



(Dieterich 1979 Ruina 1983 1, which was developed based on empirical observations and 
utilized in this study, describes the failure time of a preexisting fault (here: Beichuan 
fault system) . Furthermore, such a rate-and-state model can estimate the seismic cycle 
by taking into account a) the Mohr-Coulomb failure law, b) the slip velocity of the 
fault, and c) the history of the slip velocity: 

r — ^,<Jn, and ^ — f{v,v{t)), where (5) 

fi^ ^IO+a\n (v{t)/vo) + b In {^(t)vo/dc) (6) 

a and b are dimensionless hyper-parameters and based on empirical observations 



(Blanpied et al. 19981. /xq is the initial friction coefficient and vq is the initial slip 



velocity of the fault, ^(t) is a time dependent "state" quantity and dc is a critical slip 



distance (Dieterich 1979 Ruina 19831 



1 = 1- atHt)/dc (7) 

This model can be used to determine the clock advancement of an earthquake, for 
example, due to static triggering stress perturbations applied during a seismic cycle. 
The static triggering stress perturbation can be approximated by a heavyside function 
(eq.|8|, because the water reservoir impoundment process was abrupt in 2005 (Fig. 2). 
In detail, the clock advancement At is the period between the time of failure without 
perturbation tf and the time of failure when the static triggering load perturbed the 



background ( Gomberg et al. 1998 1 



Table 2 Model quantities for three rate-and state fault regimes describing the potential be- 
havior of the Beichuan fault system. 



model 


parameters 




rate-state model 






regime 


1 


regime 2 


regime 3 


a 




0.15 




0.005 


0.0027 


b 




0.15 




0.15 


0.15 


MO 




0.57 




0.57 


0.57 


dc m 




0.001 




0.001 


0.001 



At = tf 



\ACFS\/fi - a//iln(l - (L + 1) exp(-/ii//a)). 



(8) 



where \ACFS\ is the Coulomb failure stress of the static perturbation, normalized 
by the associated normal stress an , and L describes the load function respectively: 



exp(/ifo/a) (1 — exp{\ACFSs\/a)) — 1 if tectonic and static surface load, 
— 1 if tectonic load only. 

(9) 



Three rate-state models and three natural earthquake cycles (4kyr, 7kyr, and lOkyr) 
were chosen to estimating the clock advancement of the Wenchuan earthquake by static 
triggering stress perturbations of the artificial water reservoir behind the Zipingpu dam 
(Fig. 2). The different models and cycles will help to better understand the number 
of years by which the clock of the mainshock was advanced. Model parameters are 
summarized in the following Table [2] 
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Fig. 2. Map view of the Minjiang River Valley and the cross section of the 
double-listric fault model, including the extension of seasonal lower (L) and 
upper water level (U) in the valley (a). The construction of the Zipingpu 
dam (Z) was completed in 2004/5 | |Lei et al. 20081 1, and followed by the 
impounding of an artificial water reservoir (b). The loading time is indicated 
by to and the mainshock (black star) ruptured at tp = tf — At, while 
being clock advanced by zi t year. The water level changed seasonally with a 
wavelength of i„ = 1 year. 



4 Results and Discussion 



The period of stress perturbation was about 2.5 years, which accounts for the time 
between the start of the loading to in 2005 and the moment of failure tp in 2008. 
The 3-dimensional stress modeling results show that the static surface loading process 
during the Zipingpu water reservoir impounding altered the lithostatic stresses of the 
reverse fault conditions with a vertical minimum principal stress 0-3 and a horizontal 
maximum principal stress cti (Fig. la). In fact, static stress perturbations brought BTl 
closer to failure [ACFS > 0) and shifted BT2 and the shallow dipping root of BTl 
away from failure {ACFS < 0) as shown in Figure Ic. 

First, the water was at lower level (300 Mt) when the mainshock ruptured on May 
12, 2008. This static load induced shear stresses of >1 kPa and normal stresses > |-2| 
kPa (compression, if a <0) at <17 km depth on BTl beneath the artificial water reser- 



voir (see Section 3.1 1. Given the peak load at upper water level (1.10 Gt), shear stresses 
increased by >3 kPa and normal stresses increased by > |-6| kPa. Furthermore, the 
shear stress Art and the effective normal stress Aa'^, in turn, changed the Coulomb 
failure stress (see eq. [3| and brought BTl closer to failure in <17 km depth. Again, it 
should be emphasized that Aa'n also changed due to the poroelastic response (destabi- 
lization) of BTl to the load of the Zipingpu water reservoir (eq. 4). BTl responded, in 
particular, on steep dipping fault segments (>60°) within the low permeable Protero- 
zoic granitic rocks, as shown in Fig. 3. On the other hand, the influence of the pore 
pressure diffusion is inflnitesimal small under undrained conditions in a compressive 



tectonic regime. This fact has been previously discussed ( Klose 2008 Klose 2009 Deng 
|et al. 2010] [Gahalaut and Gahalaut 2010[ ). 



prior (correlation -0.18) 
1st cyle (correlation +0.69) 
2st cyle (correlation +0.62 i 
3st cyle (correlation -0.1 2) i 




Fig. 3. Relationship between the water level change in the Minjiang river 
valley and the earthquake event rate, recorded in the vicinity of the water 
reservoir. The correlation coefficients are provided for each period prior and 
during the surface loading process. 

Second, from October 2005 until the mainshock nucleated in May 2008, ACFS 
exceeded triggering stresses of 4 kPa at <12 km depth. This stress level is critical, 



10 



since it is daily generated by tidal elongations of the earth due to both the moon and 
the sun (see Appendix). It has been empirically shown that tidal stress changes have 
weak or random effects on triggering medium- to large-size earthquakes (Mauk and 



Kienle 1973 Klein 1976 Beeler and Lockner 2003a Cochran et al. 20041. Thus, it can 



be anticipated that any triggering stress perturbation in the earth's crust must exceed 
at least stress levels of 1-10 KPa resulting from the tidal elongation of the earth. 
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Fig. 4. Relationship between the coseismic slip s observed during the main- 
shock | |Wang et al. 2008a[ | and the clock advancement At of static triggering 
stress perturbations as result of the surface loading processes of the artificial 
water reservoir (see Fig. 2). s and At are sampled from non-interpolated finite 
elements of the double-listric Beichuan fault beneath the reservoir. Zii-values 
are based on a rate-and-state dependent friction law (see supplementary 
material). Positive/Negative Z\t- values indicate a fault de/stabilization. 
Correlation coefficients between s and At are indicated by k. 

Third, the Zipingpu reservoir only provided the trigger for the limited, initial rup- 
ture on the thrust beneath the reservoir (Fig. Ic), which grew to engage the natural 
accumulated stress Ag along the full Beichuan system (BTl-3) of about 250 km length 
(Fig. lb). Ag had been built up by natural tectonic loading during the seismic cycle 
along BTl-3 and the entire Longmen Shan. Moreover, the water reservoir also influ- 
enced the initial rupture propagation. In the first 20 seconds, the rupture propagation 
on a 32 km long and 13±3 km wide BTl-fault segment was directed upwards toward 
the Minjiang River Valley directly to the surface loading area of the Zipingpu water 
reservoir (Fig. lb). Then, the rupture propagation process changed to a right-lateral 
NNE strike-shp rupture mechanism ( [Zhang et al. 2009| |Wang et al. 2008a[ ). Thus, al- 
though expected and observed fault geometries differ, the expected seismic magnitude 
still falls, with an underestimation, into the statistical uncertainty of the observed 
magnitude. Last, modeling results based on a rate-and-state dependent friction law 
(see Section |3.2[ ) suggest that the clock advancement At oi the mainshock varies for 
different earthquake cycles and model regimes (Tab. [3|. As discussed by Gomberg et 
al. 2000, the smaller a is relatively to &, the more the fault becomes Coulomb-like. This 



11 



might be true at 11-19 km depth, where faults are generally locked in the Longmen 
Shan region (Wang et al. 2008b I. Thus, regime 3 in Table [s] might be more suitable 



to model the seismic cycle of the Beichuan faults, due to low deformation rates and 
the stiff behaving Pengguan Massif (Fig. 1). Static perturbations of the water reservoir 
more likely produced <41 years of equivalent tectonic loading in a 4 kyr earthquake 
cycle or <57 years in a 10 kyr cycle. The more (less) the surface loading of the water 
reservoir perturbed BTl and BT2 in its near vicinity (15 km radius) the more (less) 
it contributed to the coseismic slip observed on both faults during the main rupture 
(Fig. 4). This contribution, however, was <1%. Although the static triggering seems 
to contribute only up to 6 decades of equivalent tectonic loading, it strongly correlates 
with the coseismic slip observed on BTl and BT2 beneath the reservoir (Fig. 4). It 
should be noted that seasonal water level changes might advance the clock even fur- 
ther. According to the rate-and-state dependent friction law by Gomberg et al. (1998), 
dynamic stress perturbations applied late in the earthquake cycle, would advance the 
clock by several hundred years. Problematic, however, is that the friction law approxi- 
mates the dynamic load as a box function instead of a wave function, which results in 
a non-optimal solution. Thus, dynamic loads are not considered in this study. 



Table 3 Expected clock advancement times At (in years) for the 2008 Wenchuan due to 
static triggering stress perturbations A CFS of about 2 kPa at hypocentral depth and 9 kPa 
at 4 km beneath the Zipingpu water reservoir. Results are based on the rate-and-state model 
suggested by Gomberg et al. (1998) with respect to a) three earthquake cycles and b) three 
model regimes describing the potential behavior of the Beichuan fault. 
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5 Conclusion 

This study suggests that surface loading of at least 320 million metric tons of wa- 
ter, which accumulated in the Zipingpu water reservoir in the Minjiang River Valley 
between 2005 and 2008, most likely triggered advanced China's Wenchuan M7.9 earth- 
quake of May 12, 2008, while enhancing a reverse-fault rupture propagation within 
the first 20 seconds. Specifically, 3-dimensional geomechanical modeling results based 
on a rate-and-state dependent friction law show that static triggering stresses brought 
parts of the Benchuan thrust fault system nearby the water reservoir closer to failure 
and advanced the clock of the mainshock by up to six decades. Conversely, other fault 
segments, directly beneath the reservoir, were brought away from failure due to the 
weight of the reservoir. 

The estimated slip of equivalent tectonic loading that was produced by the water 
reservoir on the Beichuan fault system indicates a strong correlation with the observed 
coseismic slip during the mainshock of May 12, 2008. The highest correlation coefficient 
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of 0.93 was found for an earthquake cycle of 10 kyr. Moreover, correlations become 
weaker with decreasing recurrence time of the earthquake cycle. This confirms results 
of previous paleoseismicity studies, which show evidence that the Longmen Shan and, 
in particular, the Pengguan Massif is characterized by a >7 kyr seismic cycle for major 
M>7 earthquakes. 
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7 Appendix 



Everyday, the moon and sun cause tidal elongations on earth (Bartels 19851. These 
daily pull/push effects, in turn, induce stabilizing and destabilizing stresses on preex- 
isting fault zones in the earth's crust and are independent from any geological forces 
on earth, including endogenous forces (e.g., volcanism, tectonics) and exogenous forces 
(e.g., erosion, sedimentation). Moreover, it has been shown that tidal stress changes 



have weak effects on triggering medium- to large-size earthquakes (Beeler and Lock- 



ner 2003a I. Thus, it can be anticipated that any triggering stress perturbation on the 
earth's crust must exceed at least stress levels resulting from the tidal elongation of 
the earth. Analytically, it can be shown how high these tidal stress changes are. 

Let's assume, Fq is the earth's gravitational potential, which results from both at- 
traction force and centrifugal force of the earth and moon/sun. Tidal forces V, however, 
deform Fq: 

F ^ Fq-V (10) 

This results in a vertical surface displacement ^ with respect to the average value 
of the gravitation acceleration on earth g — 9.798 ms^ . 



Tidal forces change with geocentric zenith distance 6 from the moon/sun (Bartels 



19851, whereas the main term of the tidal potential is 



V. (I) 4 (cos 2. + 1) (12) 

with the lunar tidal constant Gi = 2.6206 m s^ and the solar tidal constant Gs 
= 1.2068 m s~ , the radius of the earth r^ and the mean radius of the earth Te = 
6371.221 km. Assuming the earth is a sphere {r^ — Y'e), the general form of vertical 
surface displacement is 

^=^(cos2e+0. (13) 
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The displacement due to the moon and sun is 

Ci = 0.267m (cos 261+ -j, ^s = 0.123m( cos26l + - j . (14) 

Thus, the ^; and ^^ are amplified at the zenith (^ = 0°) by 0.356 m and 0.164 m. 
On the other hand, they are depressed at the nadir (6 = 90°) by 0.178 m and 0.082 
m. The peak-trough difi^erence for the moon is 0.534 m and 0.246 m for the sun. 

Both vertical displacements induce maximal shear stresses r and normal stresses 
a„ on preexisting faults (e.g., dipping 45°) in the earth's crust with an average shear 
modulus of about G = 30 GPa and friction angle of, let's assume, of 29°: 

r = J- 2G, (15) 

rE 

CTn = -^, (16) 

tan0 

Thus, maximum induced stresses that need to be exceeded by any additional trig- 
gering stress (e.g., due to surface loading) are: 

by the moon r; — 5.03 kPa and (j„ ; = 9.07 kPa, 

by the sun Ts = 2.32 kPa and an,s ~ 4.18 kPa. 
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